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La presente invention concerne une methode pour detecter et controler la 
formation d' hydrates en tout point d'une conduite oii circulent des fluides p&roliers 
polyphasiques. 

Etat de la technique 

10 Produire des hydrocarbures depuis des gisements en eaux profondes (1500 - 

3000 m) souleve un grand nombre de difficultes, essentiellement a cause des 
pressions eUevees et des temperatures basses qui y regnent. Ces difficultes se 
rencontrent dans tous les domaines du savoir-faire p&rolier : le forage et les 
interventions sur puits, les installations de traitement (FPSO, FWHP, etc.) et leurs 

15 ancrages, les technologies sous-marines (pompage polyphasique, separation sous- 
marine des phases des fluides petroliers), les risers et ombilicaux ou le «flow 
assurance ». 

Dans le domaine de la production proprement dit, l'exploitant est tenu de 
garantir le bon ecoulement malgre de nombreux obstacles comme la formation de 
20 bouchons de fluide ou « slugging », les dep6ts mineraux et organiques ou bien 
encore la formation d'hydrates. En offshore profond les pressions plus elevees et les 
temperatures plus basses favorisent la formation d'hydrates. Les interruptions de 
production durant lesquelles les fluides refroidissent accroissent encore les risques de 
formation d'hydrates. 
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Une detection tardive de la presence d' hydrates notamment par 
meconnaissance des conditions de leur formation dans les conduites, peut entrainer 
des problemes de production couteux : blocage complet des conduites de production 
et/ou de transport. II est done important pour l'exploitant, de disposer de moyens 
5 d' evaluation des risques de toute nature, de fa$on a mettre en place des techniques de 
regulation et d' intervention et ainsi de pouvoir garantir un bon ecoulement des 
fluides en circulation. 

Pour proteger et isoler thermiquement les conduites sous-marines, on les 
enferme dans un tube exterieur resistant a la pression hydrostatique (« souvent 

10 designe par « bundles » par les gens de l'art). Le plus souvent, on regroupe plusieurs 
conduites pour former des faisceaux. Ces conduites regroupees ont souvent des 
fonctions differentes. Elles servent & la production d'hydrocarbures, a l'injection de 
fluides dans le gisement ou reservoir, a Tinjection de gaz dans les conduites (« gas 
lift »), k la circulation de fluide de rechauffage, etc. Dans l'espace entre elles, on 

15 interpose par exemple un calorifugeage a faible conductivity thermique laisse a la 
pression atmospherique ou mis sous vide, avec des cloisonnements places a 
intervalles reguliers pour des raisons de securite. 

Avec un outil efficace devaluation en continu des risques de formation 
dehydrates et autres depots, l'exploitant peut intervenir soit par le biais d'un 

20 rechauffage des conduites soit par l'injection d'additifs. Cependant un tel outil est 
difficile a mettre en place et couteux en temps de calcul surtout si on doit tenir 
compte de la composition detaillee des fluides petroliers. Pour etudier plus 
facilement leur comportement, il est connu de les decrire comme une combinaison 
d'un nombre de composants ou pseudo composants beaucoup plus reduit que le 

25 nombre reel de constituants. On parle alors d'une composition regroupee par 
opposition a une composition detaillee. 

Par le brevet FR 2 753 535 et la demande de brevet FR 0003867 du 
demandeur, on connait des methodes. pour predire la temperature d' apparition de 
depots tels que des cires ou paraffmes dans des bruts petroliers comportant le 
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regroupement ou « lumping » de tous leurs constituants en un nombre plus reduit de 
pseudo composants representant chacun plusieurs classes d'hydrocarbures et dont les 
parametres physico-chimiques sont determines par combinaison des parametres 
correspondants d'un certain nombre d'hydrocarbures purs regroupes dans une base 
5 de donnees. Les pseudo composants de cette formulation regroupee, sont appliques a 
un module thermodynamique permettant de determiner differents parametres 
indicatifs des conditions de formation des depots. 

Differents logiciels de simulation existent sur le marche peimettant de 
modeliser le comportement de fluides polyphasiques en circulation dans des 
10 pipelines. 

Par les brevets ou demandes de brevet suivants : US 5 550 761, FR 2.756.044 
(US 6 028 992) et FR 2 756 045 (US 5 960 187), FR 00/08 200 et FR 00/09 889 du 
demandeur, et aussi par les publications suivantes : . 

- Faille I. et Heintze E., "A rough finite volume scheme for modeling two-phase 
15 flow in a pipeline", Computers & Fluids 28 (1999), et 

- Pauchon C. et al, "TACITE: a comprehensive mechanistic model for two-phase 
flow", 6th BHRG Multiphase International Conference, Cannes, France, June 
(1993), 

on connait notamment le code de simulation TACITE (marque deposee) qui simule 
20 le comportement transitoire des fluides polyphasiques en circulation en prenant en 
compte a tout instant leur composition. Avec un tel code previsionnel, Tingenieur de 
production est en mesure de definir les caracteristiques des reseaux d'acheminement 
de fluides, notamment d'hydrocarbures : conduites, separateurs, vannes, systemes de 
controle, etc. Un des objectifs importants dans ces operations de simulation, c'est de 
25 predire avec precision les caracteristiques des ecoulements transitoires lors de 
situations telles que : variation du debit d'entree, depressurisation a la sortie, arrets 
de production et redemarrage, nettoyage des conduites (« pigging »), apparition et 
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propagation de bouchons liquides, formation de depots susceptibles de se former 
dans certaines conditions thermodynamiques, etc. 

Les modes d'ecoulement de fluides polyphasiques dans des tubes sont 
extremement varies et complexes. Les ecoulements diphasiques, par exemple, 

5 peuvent etre stratifies, la phase liquide s'ecoulant dans la partie inferieure de la 
conduite, ou intermittents avec une succession de bouchons liquides et gazeux, ou 
bien encore disperses, le liquide etant entraine sous forme de fines gouttelettes. Le 
mode d'ecoulement et ie glissement entre les phases varient notamment avec 
Tinclinaison des conduites par rapport a l'horizontale et dependent du debit, de la 

10 temperature etc. 

Le code TACITE, on le rappelle, est adapte a suivre les variations de la 
composition des melanges d'hydrocarbures. Comme deja decrit dans les brevets 
precites du demandeur, il y a une equation de conservation de masse pour chaque 
pseudo-composant. II y a aussi une equation de quantite de mouvement et une 
15 equation pour l'energie du melange. Les equations correspondantes sont rappelees ci- 
apres, apres une definition des symboles et notations: 

• V vitesse absolue de phase ; 

• U vitesse superficielle ; 

• R fraction volumetrique par phase ; 
20 • p densite de phase ; 

• H enthalpie de phase ; 

• E energie interne ; 

• P pression du melange ; 

• Tyy frottement parietal ; 

25 • Q w flux de chaleur parietal entre le pipeline et le milieu environnant ; 

• 0 inclinaison du pipeline ; 

• g gravite ; 

• S section du pipeline ; 
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x fraction massique du composant dans une phase donnee ; 
p notnbre de phases (de 1 a 3) ; 
N nombre de composants ; 
k indice de phase ; 
i indice de composant ; 
m indice de melange, 

L' equation d'equilibre de masse pour chaque composant est : 

±\£ s{p,R t A + ±\± S (p k R t x>V k )}-0 0) 
dt [t? L J o x L m J 

• L'equation d'equilibre de la quantite de mouvement du melange est 
± S{p k R k V k )Y^± s(p k R k V k 2 )yPyS(T w - Pm g sm6). (2) 

• L'equation d' 6quilibre d'6nergie du melange est : 
± S(p k E k )^+j^\t s(p k R k V k H k )yS{Q w -p m gU m sm9). (3) 



dt 



d_ 
dt 



Schema numerique 



Le schema num6rique est conservatif et non dissipatif. D fournit un bon bilan 
15 de masse et d'energie en tout point de la conduite et a tout instant. Un schema mixte 
implicite/explicite est utilis6 pour optimiser le temps de calcul et la capacite de suivi 
des fronts des ondes de taux de vide, ce qui est particulierement important dans le cas 
ou la configuration du terrain suivi par le pipeline ou sa configuration propre, 
favorise la formation de bouchons de liquide avec propagation d' ondes de vide dans 
20 les deux directions opposees du pipeline, phenomenes que les specialistes designent 
couramment par "terrain slugging" ou "severe slugging". 
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Module thermodynamique 

Le code de simulation TACITE comporte un "flash" thermodynamique 
integre i.e. un sous-programme integre de calcul des proprietes thermodynamiques 
(equilibre liquide-vapeur, composition de chacune des phases) a l'aide d'une equation 
5 d'etat. Ce flash realise des calculs d' equilibre thermodynamique di et tri-phasiques 
pour des melanges d'hydrocarbures incluant de 1'eau. Les equations d'etat cubiques 
de Peng-Robinson (1) et de Soave-Redlich-Kwong (2) sont utilisees pour modeliser 
les proprietes thermodynamiques a Tequilibre des phases. Ces equations sont 
definies dans les publications suivantes : 

10 - D.-Y. Peng D.Y. et al, "A new two-constant equation of state", Ind. Eng. Chem. 
Fund. 15, 59 (1976) ; 

- Soave G., "Equilibrium constants from a modified Redlich-Kwong equation of 
state", Chem. Eng. Sci. 27, 1197 (1972). 

Dans les deux cas, les volumes moleculaires peuvent etre corriges par la 
15 methode de Peneloux decrite dans la publication suivante : 

Peneloux A. et al, "A consistent correction for Redlich-Kwong-Soave 
volumes", Fluid Phase Equilibria 8, 7 (1982). 

Le module dit flash assure un suivi precis de la composition des fluides aussi 
bien dans Tespace que dans le temps, durant toute la simulation. 

20 Ce suivi compositionnel rend le code TACITE particulierement apte a predire 

avec precision le risque de la formation dehydrates, comme on le verra par la suite. 

Module hydrodynamique 



Le module hydrodynamique calcule le regime d'ecoulement, la vitesse de 
glissement entre phases et les termes de frottement. 
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La methode selon ^invention 

La methode selon l'invention permet de detecter en continu en tout point 
d'une conduite oil circule un melange polyphasique de fluides petroliers, les 
conditions thermodynamiques de formation d'hydrates. Elle comporte l'utilisation 

5 d'un module hydrodynamique mecanistique et d'un module thermodynamique 
compositionnel integre pour definir les proprietes des phases et la resolution 
d'equations de conservation de la masse, de conservation de la quantite de 
mouvement et de transfert d'energie au sein du melange, en considerant que le 
melange de fluides est sensiblement a Tequilibre a chaque instant, que la composition 

10 du melange polyphasique est variable tout le long de la conduite et que la masse de 
chaque constituant du melange est definie globalement par une equation de 
conservation de la masse sans consideration de son etat de phase. Les fluides 
petroliers sont regroupes en un nombre limite de pseudo-composants regroupant 
chacun certaines fractions d'hydrocarbures. 

15 La methode est caracterisee en ce que Ton detecte les conditions 

thermodynamiques de formation d'hydrates apres avoir opere un regroupement 
particulier des fluides petroliers en pseudo-composants de fa$on a isoler les 
composants formateurs d'hydrates, avec definition pour chacun d'eux d'une fraction 
massique et d'un certain nombre de grandeurs physiques caracteristiques, et on 

20 applique aux modules thermodynamique et hydrodynamique, des donnees relatives k 
ces pseudo-composants selectionnes de maniere a determiner en tout point la 
temperature (Td) de dissociation des hydrates. 

L'invention concerne aussi une methode pour controler en continu et en tout 
point d'une conduite ou circulent des fluides petroliers polyphasiques, la formation 
25 d'hydrates. Elle comporte la detection des conditions thermodynamiques de 
formation des hydrates comme defini ci-dessus, l'utilisation d'un organe de controle 
pour comparer la temperature des fluides petroliers a cette temperature de 



8 



dissociation (Td) et T application commandee par cet organe de controle, de mesures 
destinees a combattre la formation des hydrates. 

Suivant un premier mode de mise en oeuvre, on utilise des moyens de 
chauffage associes a la conduite pour elever la temperature des fluides petroliers au- 
5 dessus de la temperature de dissociation. Dans le cas ou la conduite est incluse avec 
au moins un deuxieme conduit dans un tube isole du milieu exterieur, on utilise le 
deuxieme conduit pour y faire circuler un fluide chaud. 

Suivant un deuxieme mode de mise en oeuvre, sur commande de Torgane de 
controle (C), on injecte des inhibiteurs d'hydrates dans la conduite. 

10 

Presentation sommaire des figures 

D'autres caracteristiques et avantages de la methode selon Tinvention, 
apparaitront a la lecture de la description ci-apres d'un exemple non limitatif de 
realisation, en se referant aux dessins annexes ou : 

15 - la Fig.l montre un organigramme compare des operations de simulation dans le 
cadre de Tapplication du code TACITE selon que Ton etudie ou non le risque de 
formation d' hydrates ; 

- la Fig.2 montre un organigramme de deroulement des operations de simulation 
qui se repetent en tous les points d'un maillage de pipeline ; 

20 - la Fig.3 montre schematiquement un mode de controle que Ton peut appliquer 
pour mettre en oeuvre toute intervention visant a combattre la formation 
d' hydrates ; 

- la Fig.4 montre la section transversale d'un tube le long duquel courent plusieurs 
ensembles de tubes ou lignes de circulation de fluides ; 
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la Fig.5 montre une comparison des temperatures de dissociation d'hydrates 
calculees aux temperatures exp&imentales dans le cas du methane. 
L'histogramme est etabli a partir de 158 points experimentaux ; 

la Fig.6 montre une comparaison des temperatures de dissociation d'hydrates 
calculees aux temperatures experimentales dans le cas de l'ethane. 
L'histogramme est etabli a partir de 126 points experimentaux ; 

la Fig.7 montre une comparaison des temperatures de dissociation d'hydrates 
calculees aux temperatures experimentales dans le cas du propane. 
L'histogramme est etabli a partir de 93 points experimentaux ; 

la Fig.8 montre une comparaison des temperatures de dissociation d'hydrates 
calculees aux temperatures experimentales dans le cas de l'isobutane. 
L'histogramme est etabli a partir de 43 points experimentaux ; 

la Fig.9 montre une comparaison des temperatures de dissociation d'hydrates 
calculees aux temperatures experimentales dans le cas du dioxyde de carbone. 
L'histogramme est etabli a partir de 102 points experimentaux ; 

la Fig. 10 montre une comparaison des temperatures de dissociation d'hydrates 
calculees aux temperatures experimentales dans le cas de T azote. L'histogramme 
est etabli & partir de 57 points experimentaux ; 

la Fig. 11 montre une comparaison des temperatures de dissociation d'hydrates 
calculees aux temperatures experimentales dans le cas du sulfure d'hydrogene. 
L'histogramme est etabli a partir de 64 points experimentaux ; 

la Fig. 12 montre une comparaison des temperatures de dissociation d'hydrates 
calculees aux temperatures experimentales dans le cas de melanges 
d'hydrocarbures. L'histogramme est etabli a partir de 55 points experimentaux 
provenant de 8 melanges differents ; 
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- la Fig. 13 montre schematiquement un exemple de raccordement d'un pipeline 
avec un riser avec un point d' injection de gaz (pour gas lift) ; 

- la Fig. 14 montre dans un diagramme bidimensionnel longueur de conduite- 
temps, un exemple de domaine que l'on peut calculer par la methode d'apparition 

5 d'hydrates dans Y agencement de la Fig. 13 ;et 

- la Fig. 15 montre dans un diagramme bidimensionnel pression-temperature deux 
courbes locales d'apparition d'hydrates dans les conduites, selon que Ton tient 
compte ou non d'une modification de la composition du melange en circulation. 

10 DESCRIPTION DETAILLEE 

Presentation du module thermodvnamique 

Structure des hydrates 

Les hydrates sont des composes solides constitues de molecules 
d'hydrocarbures (methane, ethane, dioxyde de carbone etc.) piegees a haute pression 
15 et basse temperature a 1'interieur de cages approximativement spheriques formees 
par un reseau cristallin tridimensionnel de molecules d'eau. Cette structure est 
metastable en l'absence d'hydrocarbures et necessite un taux de remplissage minimal 
pour assurer sa stabilisation. A la difference du reseau cristallin de la glace, celui des 
hydrates peut demeurer stable au-dessus de 0°C. 

20 Des mesures par diffraction aux rayons X realisees dans les annees 50 ont 

permis de mettre en evidence deux structures cristallines des hydrates selon la nature 
des molecules d' insertion presentes. Dans les deux cas, il s'agit de structures m 
cubiques organisees en reseaux de petites et grandes cavites. Le tableau ci-apres 
montre T occupation possible des differentes cavites selon la nature des molecules 

25 formatrices d'hydrates mises en jeu (+ : occupation possible de la cavite par la 
molecule ; - : taille de la cavite insuffisante pour contenir la molecule). 
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Le reseau de la structure I est de type cubique face centree. Chaque maille 
contient 46 molecules d'eau organisees en : 

2 petites cavites de 3,95 A de rayon, 
6 grandes cavites de 4,33 A de rayon. 

Le cristal de la structure II a le reseau du diamant. Une maille contient 136 
molecules d'eau et 24 cavites : 

16 petites cavites de rayon 3,91 A, 
8 grandes cavites de rayon 4,73 A. 

Une troisieme structure, la structure H, a ete egalement identifier Son reseau 
cristallin, de type hexagonal, est constitue de 34 molecules d'eau et de 6 cavites : 

2 petites cavites de 3,91 A de rayon, 

3 cavites intermediaires d'un rayon de 4,06 A, 
1 grande cavite de 5,71 A de rayon. 
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Cette structure est beaucoup plus rare que les deux precedentes, elle se 
rencontre principalement dans le cas d'alcanes ou d'alcenes cycliques. 

ThermodvnamiQue de Tequilibre de phases en presen ce d'hvdrates 

Tout calcul de prediction d'hydrates est base sur le critere d'equilibre de 1'eau 
5 en phase hydrate et en phase aqueuse, c'est-a-dire sur l'6galite du potentiel chimique 
de l'eau dans la phase hydrate et du potentiel chimique de l'eau dans la phase 
aqueuse : 

f,^[p,T) = ^{P,T). (1) 

Selon la theorie statistique de Van Der Waals et Platteeuw decrite dans la 
10 publication suivante :Van der Waals, J.H., et al., "Clathrate Solutions", Adv. Chem. 
Phys., 2, 1, (1959), 

le potentiel chimique de l'eau dans la phase hydrate peut s'ecrire : 
^{P.T^^-^P^ + RTZv, Infl-ZflA (2) 



avec : 



15 - /J^r tc " vide : potentiel chimique de l'eau dans le reseau cristallin en Tabsence de 
molecule d'hydrocarbure, 

- v. : rapport du nombre de cavites de type i (petite ou grande) sur le nombre de 
molecules d'eau dans un monocristal, 

- : probabilite pour qu'une cavite de type i soit occupee par une molecule 
20 hydrocarbonee de type k. D'apres la theorie de 1'adsorption de Langmuir, cette 

probabilite peut se calculer par : 
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Ctkfk 



ou fk est la fugacite du composant k et C» est le coefficient d' adsorption de 
Langmuir d'une molecule k dans une cavite de type i. 

Les coefficients de Langmuir peuvent etre exprimes a partir du potentiel 
d'interaction entre les molecules d'eau et les molecules de gaz incluses dans les 
cavit6s de la structure cristalline d'hydrate : 



oil k est la constante de Boltzmann et W(r) l'energie potentielle d'interaction entre les 
molecules d'eau formant la cavit<5 et la molecule de gaz se trouvant a une distance r 
du centre de la cavit6. II existe plusieurs formes de potentiel, les plus utilisees sont le 
potentiel de Lennard- Jones, le potentiel de Kihara et le potentiel a puits carr6. 

Le potentiel de Lennard- Jones entre deux molecules s'ecrit : 





avec : 



e : profondeur du puits, 

r : distance entre les centres des molecules, 

o : diametre d'exclusion, r tel que W( a) = 0. 



Le potentiel de Kihara est identique a celui de Lennard-Jones, mais les 
molecules sont considerdes comme des spheres de rayon non nul. Son expression 
analytique est la suivante : 



iui uupui 



14 



W(r) = 



si r< 2a 



(6) 



W(r) = 4e 




si r > 2a 



avec : 



e : profondeur du puits, 

a : tel que W( a+2a) = 0, 

a : rayon moyen des 2 molecules, 

r : distance entre les centres des molecules. 

Enfin, le potentiel a puits carre est donne par : 

W(x) = oo si 0<x<<7 

W(x) = -6 si a < x < ao (7) 
W(x) = 0 si x>ao 

Les differents parametres figurant dans les expressions de ces potentiels sont 
ajustes sur des donnees experimentales d'equilibre d'hydrates. 

Le potentiel chimique de l'eau dans la phase aqueuse s'ecrit : 



ou y est le coefficient d'activite de l'eau, x la fraction molaire de l'eau dans la phase 
aqueuse et ;U° u (P,r) le potentiel chimique de l'eau pure. C'est dans ce calcul de 

l'activite de l'eau en phase aqueuse que sont entre autres pris en compte les effets des 
sels et des alcools. 

En egalisant les equations (2) et (8), la condition d'equilibre recherchee (1) 

s'ecrit : 



/C ueuse (^) = ^(P,T)+RT\n{ycl 



(8) 



nL(P,T)-fi, 



hydrate -vide 
'eau 




10 
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En fixant la pression, cette equation nous pennet d'acceder a la temperature 
d'equilibre ou inversement, en fixant la temperature nous obtenons la pression 
d'equilibre. 

Les differents modules de prediction d'hydrates sont tous bases sur le 
formalisme decrit ci-dessus. Ces modules vont essentiellement se differencier dans le 
choix du potentiel d'interaction eau - hydrocarbure, dans le choix de 1'equation 
d'etat utilisee pour le calcul des fugacites ainsi que dans le module de calcul de 
Tactivite de l'eau en phase liquide. Dans ce qui suit, nous presentons plus 
precisement le module que nous avons adapte pour les besoins de TACITE. 

Module implemente dans TACITE 

Le module de prediction d'hydrates implemente dans TACITE s'inspire du 
module de Munck tel que decrit par Munck, J. et al, "Computations of the formation 
of gas hydrates", Chem. Eng. Sci 43, 2661 (1988). II permet de calculer la 
temperature de dissociation des hydrates h pression donnee dans le cas d'un gaz, 
d'un liquide ou d'un melange gaz - liquide diphasique en contact avec de Teau 
liquide. 

Constantes de Langmuir : 

Ce module utilise un potentiel W(r) de type puits carre, ce qui conduit a une 
expression de la constante de Langmuir du compose k dans une cavite de type i sous 



la forme : 




Aik ,Bik. 
— exp( — ) . 



(10) 
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Les valeurs des parametres A et B des differents composes formant des hydrates sont 



donnees dans le tableau ci-apres : 







Petite cavite 


Grande 


cavite 


Molecule 


Structure 


A v 10 3 


B 


Ax 10 3 


B 






^XW ami ) 


(X) 


(K/atm) 


(K) 


N 2 


I 


1 617 


2905 


6,078 


2431 




II 


0 1742 


3082 


18,00 


1728 


C0 2 


I 


0 2474 


3410 


42,46 


2813 




n 


0 0845 


3615 


851,0 


2025 ! 


H 2 S 


1 


0,0250 


4568 


16,34 


3737 




n 


0,0298 


4878 


87,2 


2633 


c, 


1 


0,7228 


3187 


23,35 


2653 




n 


0,2207 


3453 


100,0 


1916 


c 2 


1 


0,0 


0,0 


3,039 


3861 




11 


0,0 


0,0 


240,0 


2967 


c 3 


11 


0,0 


0,0 


5,455 


4638 


iC 4 


11 


0,0 


0,0 


189,3 


3800 


nC 4 


11 


0,0 


0,0 


30,51 


3699 



Calcul des fugacites : 

Les fugacites des differents constituants presents dans le melange, formateurs 
d'hydrates ou non, sont calculees a Taide de l'equation d'etat de Soave-Redlich- 
Kwong deja citee. Un tel calcul se fait par l'appel au module de flash a pression et 
temperature fixees integr6 dans TACITE. 

Activite de Teau en phase aaueuse : 

Le module de Munck neglige la solubilite dans la phase aqueuse des 
hydrocarbures, de l'azote et du sulfure d'hydrogene. En revanche, la solubilite du 
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dioxyde de carbone dans l'eau ne peut etre negligee et est calculee a l'aide de la 
constante de Henry du C0 2 dans l'eau : 



x co 2 ~ 



fco 2 
H co 2 ' 



(11) 



°™ x co z est la fraction molaire de C0 2 dans la phase aqueuse, f co ^ la fugacite du C0 2 
et H COz la constante de Henry du C0 2 dans l'eau donnee par : 



lnH COj =H A +£?- + H c lnT + H D T. 



(12) 



Les valeurs des coefficients utilises pour le calcul de la constante de Henry du 
C0 2 figurent dans le tableau ci-apres. 



(atm) 


H B 
(atm.K) 


H c 
(atm/InK) 


H D 
(atm-K" 1 ) 


160,27 


-8764,5 


-21,726 


1.1055.10" 4 



Constantes de reference : 

La difference de potentiel chimique entre l'eau liquide pure et l'eau dans la 
phase hydrate vide (premier membre de l'equation 9), peut s'ecrire, moyennant 
quelques approximations : 



J RT 2 }r? • 



RT 



RT n 



(13) 



Dans cette equation, A|t 0 represente la difference de potentiel chimique de 
l'eau dans le reseau hydrate vide et dans l'eau liquide a T 0 = 273,1 5 K. AH 0 est la 
difference d'enthalpie correspondante, AC p la difference de capacite calorifique et 
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m 



10 



AV la difference de volume. La pression P 0 est la pression de vapeur a r o ; etant 
donne que P 0 est tres petit compare a P , on utilise P 0 = 0 . 

La dependance en temp6rature du terme PV est moyennee par : 
r = (r + 273,15)/2. (14) 

et A# 0 sont tous deux considered comme etant ind6pendants de la 
pression. AV est considere aussi comme independant de la temperature, alors que la 
dependance en temperature du terme enthalpique est prise en compte par le biais 
d'une difference de capacite calorifique constante AC p . Les differentes constantes 
figurant dans liquation (13) sont presentees dans le tableau suivant : 



Proprietes 


Unites 


Structure I 


Structure II 




J.mol' 1 


1264 


883 




Lmol" 1 


-4858 


-5201 


AV 


cm'-mol* 1 


1151 


808 


AC P 


J-mol 1 .K- 1 


39,16 


39,16 
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La combinaison des equations (9) et (13) donne l'equation finale (15) : 

^_|AH 0 +Ac(r-rJ r i (15) 

Dans cette equation les inconnues sont la pression et la temperature. Nous 
avons arbitrairement choisi de fixer la pression et de calculer la temperature 
d'equilibre des hydrates correspondante. La recherche de cette temperature, 
c'est-a-dire la resolution de l'equation (15), s'effectue de maniere iterative par une 
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m£thode de Newton-Raphson partant d'une temperature initiale arbitrairement fixee 
a 0°C. Nous avons verifie que la valeur de cette initialisation n'influe pas sur le 
resultat obtenu. 

Validation du module 

Nous avons valide le programme sur un ensemble de 700 points d'equilibre 
experimentaux issus de 36 references bibliographiques differentes. Les resultats sont 
presentes sous forme d'histogrammes ou figurent en abscisse, l'ecart entre la 
temperature de dissociation calculee et celle mesuree, et en ordonnee, le pourcentage 
de points pour lesquels ces ecarts de temperature sont compris entre les limites 
indiquees en abscisse. Les points experimentaux utilises pour realiser ces 
histogrammes correspondent a diff€rents couples pression - temperature. 

Sur les figures 5 & 11, nous comparons ainsi les resultats du calcul avec les 
donnees experimentales dans le cas de divers corps purs : methane, ethane, propane, 
isobutane, dioxyde de carbone, azote et sulfure d'hydrogene. 

15 Nous constatons que sur Tensemble des corps purs etudies, le module utilise 

donne d'excellentes predictions des conditions d'equilibre des hydrates. Pour 99% 
des points, la temperature est predite a ±1°C. Dans tous les cas, Tecart maximal entre 
les temperatures de dissociation calculees et mesurees ne depasse pas 1,6°C. 

L'histogramme concernant les donnees de melanges d'hydrocarbures, gaz 
20 naturels ou huiles, est presente sur la Fig. 12. Les donnees experimentales utilisees 
pour calculer cet histogramme proviennent de huit melanges differents, six de type 
gaz naturels et deux huiles. On observe a nouveau une tres bonne description du 
comportement de ces systemes par le module de prediction utilise. 



5 
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Simulation globale avec ou sans calcul de risque tf ap parition d'hvdrates 

L'impact se situe au niveau du regroupement des composants de base 
(lumping). Dans le cas standard, on cherche a en optimiser le nombre de fa<?on a 
reduire le temps de calcul tout en gardant une bonne representation des proprietes 
5 d'equilibre et des proprietes des phases. Dans le cas ou Ton cherche a predire les 
hydrates, il faut pouvoir suivre la composition des composants appeles formateurs 
d'hydrates et done, ces composants ne doivent pas etre regroupes avec d'autres. 

Outre les donnees se situant dans les fichiers de type « .PVT », « .MAS » et 
« TYP » decrits precisement ci-apres, les differentes donnees d' entree de TACITE 
10 sont regroupees dans 6 fichiers (cf. Fig. 1) : 

« .SCE » : description du scenario, e'est a dire de revolution dans le temps de la 
valeur des differentes conditions aux limites (debit d' entree, pression de sortie, 
ouverture de vanne,. , .) ; 

- « .TOP » : description de la topographie de la conduite (elevation de la conduite 
15 en fonction de la longueur) ; 

« .GEO » : description de la geometrie de la conduite (diametre, longueur, 
porosite) ; 

« .THE » : description des isolations thermiques de la conduite ; 
« .MSH » : description du maillage ; 

20 - « .STO » : description du type de fichier de resultats et de la frequence de 
stockage des resultats. 

Fichier «.PVT» 



Le fichier «.PVT» contient autant de lignes que de donnees 
thermodynamiques du fluide a indiquer, respectivement : 
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- un entier qui definit le type d'6quation d'etat utilisee : 3 pour Peng-Robinson ou 
4 pour Soave-Redlich-Kwong. 

un entier qui definit le nombre de phases ; il est mis a 1 et c'est le flash qui 
determine leur nombre. 
5 - le nombre de compos ants. 

- Boucle sur les composants avec indication, pour chacun d'eux, de leur : 

- masse molaire [kg/mol], 

- temperature critique [K]. 

- pression critique [Pa]. 

10 - volume critique [m 3 /mol]. 

- facteur acentrique [-]. 

- parachor [((N/m) 1/4 m 3 /mol)]. 

- volume molaire standard (inutilise). 

- facteur de correction volumique [m 3 /mol], 

15 - coefficients d'enthalpie ideale (tableau de 7 coefficients decrits ci-dessous). 

- Fin de la boucle. 

- la matrice des coefficients d'interaction binaire [-]. 

♦ 

- la temperature de calage 1 [K], Temperature de calage 2 [K]. 

- la viscosite cinematique a la temperature 1 [m 2 /s], Viscosite cinematique a la 
20 temperature 2 [m 2 /s]. 

Fichier «.MAS» 

Ce fichier comporte egalement un certain nombre de lignes pour indiquer : 

- le nombre de composants initiaux. 

- la liste du numero du groupe dans lequel le composant initial correspondant a ete 
25 mis. L'ordre est celui dans lequel se trouvaient les composants dans le fichier 

«.INP ». 

- le nombre de pseudo-composants. 

les differentes fractions massiques des composants, ordonnees comme dans le 
fichier «.PVT». 
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Fichier «/TYP» 

Ce fichier contient autant de Hgnes que de pseudo-composants definis dans le 
fichier «.PVT» et les informations sur les composants tiennent compte de Tordre 
dans lequel les composants sont ranges dans le fichier «.PVT». 

La i 6me ligne concerne Tidentificateur de la nature du i eme composant. 

La convention suivante est retenue : 

Si le composant est ni 1'eau, ni un formateur d'hydrates, son type est mis a 0. 

Si le composant est le sulfure d'hydrogene, son type est mis a 1. 

Si le composant est 1* azote, son type est mis a 2. 

Si le composant est le dioxyde de carbone, son type est mis 3. 

Si le composant est le methane, son type est mis a 4. 

Si le composant est l'ethane, son type est mis a 5. 

Si le composant est le propane, son type est mis a 6. 

Si le composant est Fiso-butane, son type est mis a 7. 

Si le composant est le n-butane, son type est mis a 8. 

Si le composant est Teau, son type est mis a 9. 

Le code de simulation TACITE utilise une methode de resolution numerique 
de type a volumes finis pour simuler les ecoulements dans les installations de 
production, les conduites etant discretisees par exemple selon la methode de maillage 
decrite dans la demande de brevet FR EN 00/08200 du demandeur. 

A tout instant, en tout point du maillage, les calculs de modelisation sont conduits 
selon Forganigramme de la Fig.2. 

S Le schema numerique est alimente par la donnee des termes de flux et 
les termes source, la modification des conditions aux limites (scenario). 

S Le schema numerique fournit les masses des constituants, la quantite 
de mouvement du melange, 1'energie interne du melange (non-isotherme). 
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' La connaissance des masses et de l'energie fournit, grace au flash 
thermodynamique, la pression, la temperature et la composition massique des phases 
a l'6quilibre. 



' Connaissant la pression et la composition massique, il est possible, 
5 1'utilisateur le souhaite (voir partie entouree d'un trait en pointilles), de : 



si 



calculer precisement grace au module de prediction d'hydrates, la 
temperature Td de dissociation des hydrates (temperature au-dessous de laquell 
des hydrates se sont formes, ils ne peuvent se dissocier) ; 



e si 



d'utiliser un controleur C de type PID (Fig.3) pour comparer la 
0 temperature calculee du fluide (Tfluide) a cette temperature Td de dissociation ; et 

d'agir si cette temperature de fluide est inferieure a la temperature de 
dissociation majoree d'un petit ecart (pour anticiper le probleme eventuel) sur Tun 
des moyens de contrdle prevu . 

' Independamment de ce contrdle eventuel, le calcul continue, la 
connaissance de 1'equilibre du melange et des proprietes des phases, associeVa la 
connaissance de la quantite de mouvement du melange permet, grace au module 
hydrodynamique, de fournir le glissement entre les phases et les frottements 
parietaux et interfaciaux. 

' II est alors possible de calculer les termes sources et les flux pour le 
prochain calcul. 

Controle en lip n^ 

Les operations de lutte contre la formation d'hydrates sont generalement tres 
couteuses. La moddisation fine permise par la m6thode selon 1'invention, permet a 
l'exploitant de connaitre a tout instant les conditions thermodynamiques regnant tout 
le long du pipeline et done de determiner les endroits ou les hydrates sont 
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suscpdb.es de se former. .. est en mesure quand >e besoin s e„ fart senur « 
seu.ement a ces moments-.,, de tnettre en omvre des actions poor combame .es 
a^ts dmydra.es. Les necessi.es dMmervemion son. p.us particuherement crmqu s 
„ cas « momemane de ,a production ou quand on precede a des injecuons de 
gaz dans le melange en circulation (aux emplacements de . gas lift ») par exemple. 

Un premier mode d'action consiste a rechauffer les fluidea en circulation de 
m aniere a les amener au-dessus de la .emperarure de dissociation dehydrates que . on 
. predire, au moins aux endroto de la conduite ou le risque es. apparu. Us condones 
pLlieres pour la production marine par grands funds son. generalemen. mc uses 
, avec Cannes hgnes dans une conduite meraUique T revetue exterieuremen, d une 
couche iso.ame. Cer ensembie de Ugnes (deaigoe par <« bund.es>> par les gens de 
,.») est reparti souven. en p.usieurs faisceaux .els que Fl. F2 (cf. Fig.4). Cheque 
faisceau comporte generalemen, une ou plusieurs lignes de transport 1, 2, pour 
fluides monophasiques ou po.yphaa.ques qui peuven, etre rassemb.ees dans un ,u be 
15 commun 4 evenrue.lemen, reconvert lui-meme exterieuremen, d'un. couche .solanre. 
L-inrerieur de cheque mbe 4 est rempli le plus souven, d'un fl uide monophas.que 
L >es P ace interieur de ,a condui.e T, entte .es different* faisceaux Fl, F2, est remph 
par exemp.e d'une mousse mermiquement iso.ante 5. L'une des Hgnes de transport 2, 
3 pent tee u«i,isee par exemp.e pour faire circuler si besom es,, un flume de 
20 rSchauffage. 

Pour modeliser les tansferts mermiques entre les condui.es petrolitas de 
Pensemb.e de hgnes e, .eur environnemen. (autres Ugnes. milieu exterieur) en tenant 
comp<e des moyens d'iso.arion thermique iniercales, on pen. se reporter par exemp.e 
a la publication suivante : 
25 Duret E e. al, "Pipeline bundles model imp.emen«ed into a multiphase flow 
moder, SPE 62948, SPE Annual Technical Conference and Exhibition, Dallas, 
Texas, 1-4 October (2000). 
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Un deuxieme mode d' action consiste a injecter dans les fluides des additifs 
inhibiteurs de type connu s'opposant a la formation d'hydrates. 

II s'agit d'un systeme de controle base, non pas sur une mesure fournie par un 
capteur reel, mais sur une mesure fournie par un capteur dit « logiciel ». 

5 Interet d'un code compositionnel pour un calcul d' hydrates 

L'utilisation d'un code compositionnel tel que le code TACITE permet en 
tout point du pipeline, un calcul exact, a la pression calculee du fluide et pour la 
composition locale calculee en cours de simulation, de la temperature de dissociation 
des hydrates, et par consequent une bien meilleure simulation des risques 

10 d'apparition d'hydrates que des codes ou la composition est supposee fixe, 
correspondant a celle que Ton injecte a l'entree de la conduite. On a pu verifier par 
exemple que l'injection de gaz au pied d'un « riser » (technique bien connue des 
gens de Tart), modifiait les conditions thermodynamiques locales d'un melange 
donne de fluides petroliers, au point que Ton pouvait predire la formation d'hydrates, 

15 alors qu'en supposant une composition fixe du melange, ce risque n'etait 
aucunement detect^. Nous illustrons ce point avec la simulation effectuee d'une 
injection de gaz au pied d'un riser se raccordant a un pipeline tel que celui 
schematise sur la Fig. 13, alimente par un melange dont la composition est donnee 
dans le tableau ci-dessous. 

20 



Components 


Composition (% mol.) 


c, 


10.53 


c 2 


2.64 


c 3 


2.63 


IC 


2.63 


NC 4 


2.63 


ic 5 


2.63 


NC 5 


2.63 


NC,o 


13.16 


NCo 


7.90 


H 2 0 


52.62 



lei uipui 
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La simulation est effectuee en considerant que le melange initial a ete modifie 
par une injection en bas du riser, de methane favorisant la formationd' hydrates. 

Dans un diagramme longueur de conduite-temps, le domaine d' apparition 
d'hydrates que Ton obtient est schematise a la Fig. 14. 

5 On a calcule et represents sur le diagramme bidimensionnel pression-temps 

de la Fig. 15, deux courbes locales de formation d'hydrates, Tune (en pointilles) en 
considerant que la composition n'a pas ete alteree par l'injection de methane, comme 
le ferait un code de simulation non compositionnel, 1' autre (en traits pleins) en tenant 
compte de la modification de composition intervenue. 

10 On verifie aisement sur la Fig. 15, que le point qui y est figure representant le 

couple (pression, temperature) a l'etat final serait considere a tort comme sans risque 
puisqu'& droite de la courbe d'hydrates calculee a la composition initiale. Le 
domaine de risque d' apparition d'hydrates, presente sur la Fig.14, serait egalement 
modifie en raison des variations de composition. 

15 
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REVENDICATIONS 

1) Methode pour detecter en continu en tout point d'une conduite ou circule 
un melange polyphasique de fluides petroliers, les conditions thermodynamiques de 
formation d'hydrates dans laquelle on utilise un module hydrodynamique 
5 mecanistique et un module thermodynamique compositionnel integre pour definir les 
proprietes des phases et l'on applique des equations de conservation de la masse et 
de la quantite de mouvement, et de transfer! d'energie au sein du melange, en 
considerant que le melange de fluides est sensiblement k 1'equilibre a chaque instant, 
que la composition du melange polyphasique est variable tout le long de la conduite 
10 et que la masse de chaque constituant du melange est definie globalement par une 
equation de conservation de la masse sans consideration de son etat de phase, et on 
regroupe les fluides petroliers en un nombre limite de pseudo-composants, 
caracterisee en ce que Ton detecte les conditions thermodynamiques de formation 
d'hydrates en : 

15 - operant un regroupement particulier des fluides petroliers en pseudo-composants 
de fa?on h isoler les composants formateurs d'hydrates avec definition . pour 
chacun d'eux d'une fraction massique et d'un certain nombre de grandeurs 
physiques caracteristiques ; et 

- en appliquant aux dits modules, les donnees relatives a ces fractions particulieres, 
20 de maniere a determiner en tout point la temperature (Td) de dissociation des 

hydrates. 

2) Methode pour controler en continu la formation d'hydrates, et en tout point 
d'une conduite ou circule un melange de fluides petroliers polyphasiques, dans 
laquelle on utilise un module hydrodynamique mecanistique et un module 
25 thermodynamique compositionnel integre pour definir les proprietes des phases et 
Ton applique des equations de conservation de la masse et de la quantite de 
mouvement, et de transfert d'energie au sein du melange, en considerant que le 
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melange de fluides est sensiblement a l'equilibre a chaque instant, la composition du 
melange polyphasique est variable tout le long de la conduite et la masse de chaque 
constituant du melange est definie globalement par une equation de conservation de 
la masse sans consideration de son etat de phase, et on regroupe les fluides petroliers 
en un nombre limite de pseudo-composants, caracterisee en ce que : 

a) l'on detecte les conditions thermodynamiques de formation d' hydrates en : 

- operant un regroupement particulier des fluides petroliers en pseudo-composants 
de facon a isoler les composants formateurs d'hydrates, avec definition pour 
chacun d'eux d'une fraction massique et d'un certain nombre de grandeurs 
physiques caracteristiques ; et 

- en appliquant aux dits modules, les donnees relatives a ces fractions particulieres, 
de maniere a determiner en tout point la temperature (Td) de dissociation des 
hydrates ; 

b) on utilise d'un organe de controle (C) pour comparer la temperature des 
fluides p6troliers a cette temperature de dissociation (Td) ; et 

c) on applique sous le controle de cet organe de controle (C), des mesures 
destinees a combattre la formation des hydrates. 

3) Methode selon la revendication 2, caracterisee en ce que Ton utilise des 
moyens de chauffage associes a la conduite (1) pour elever la temperature des fluides 
petroliers au-dessus de la temperature de dissociation. 

4) Methode selon la revendication 3, caracterisee en ce que, la conduite (1) 
6tant incluse avec au moins un deuxieme conduit (2, 3) dans un tube (T) isole du 
milieu exterieur, on utilise le deuxieme conduit (2, 3) pour y faire circuler un fluide 
chaud. 

5) Methode selon la revendication 2, caracterisee en ce que, sur commande de 
l'organe de controle (C), on injecte des inhibiteurs d'hydrates dans la conduite (1). 
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